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Abstract 

A high signal-to-noise spectrum covering the largest number of hydrogen lines (90 H2 lines 
and 6 HD lines) in a high redshift object was analyzed from an observation along the sight-line 
to the bright quasar source J2123— 005 with the UVES spectrograph on the ESO Very Large 
Telescope (Paranal, Chile). This delivers a constraint on a possible variation of the proton-to- 
electron mass ratio of A/i/^u = (8.5 it 3.6stat i 2.2syst) x 10^® at redshift Zabs = 2.059, which agrees 
well with a recently published result on the same system observed at the Keck telescope yielding 
A/x//x = (5.6ib 5.5stat =t 2.9syst) x 10~^. Both analyses used the same robust absorption line fitting 
procedures with detailed consideration of systematic errors. 

PACS numbers: 06.20.Jr, 95.85.Mt, 98.80.Es, 33.20.-t 



Experimental searches for temporal variation of fundamental constants on a cosmological 
time scale have resurfaced on the agenda of contemporary physics since the ground-breaking 
study of a possible variation of the fine-structure constant a from quasar observations, now a 
decade ago [1]. While a = e'^/Aireofic is a fundamental constant representing the strength of 
the electromagnetic force, the proton-to-electron mass ratio fi = mp/rrie is another important 
dimensionless constant to be tested for possible temporal drifts. Since the quark masses 
contribute only marginally to the proton mass (~ 10%), rrip is nearly proportional to the 
strength of the nuclear force, Aqcd- Hence, /i scales with Aqcd/'^'h, vh being the Higgs 
vacuum expectation value. Therefore, A/i//i probes the cosmological evolution of the nuclear 
vs. the electroweak sector in the Standard Model. Interestingly, various models of Grand 
Unification predict relations of the form 



where the proportionality constant is normally large, i.e. |-R^q,| ~ 10 — 40 [2], although 
model dependent [3]. Measuring any possible changes in a and fj, therefore allows a test of 
'beyond Standard' theories. 

At cosmological redshifts z > 2 molecular hydrogen is the test system for probing a drift 
in /i [4ii6j, where the time arrow is defined by A/i = /i^ — /^o- Wavelengths of spectral 
lines in the Lyman and Werner bands of H2 observed in distant objects, may be compared 
to the wavelengths A° for the same lines measured in the laboratory at high accuracy [6H8] 
via 



with Zabs the redshift of the absorbing gas cloud in the line-of-sight of the quasar and Ki 
the so-called sensitivity coefficient, expressing the sensitivity for each spectral line to a mass 
variation Afi; the Ki can be calculated to an accuracy of 1% [H El [7]. For the deuterated 
isotopologue HD, accurate laboratory data and K-coefficients are available |9]. 

Here we report on an investigation of the H2 absorption system at Zabs = 2.06 toward the 
quasar J2123— 005 making use of the Ultraviolet and Visible Echelle Spectrograph (UVES) 
at the Very Large Telescope (ESO, Paranal, Chile). The same object was recently inves- 
tigated with the HIRES instrument at the Keck Telescope [TO], so that now, for the first 
time, an accurate /i-variation analysis becomes available from two independent telescopes. 
Moreover, J2123— 005 is the brightest quasar with H2 absorption, with R = 15.8 mag. An 



A/i//i = R^a^a/a 
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FIG. 1: (Color online) Spectrum of J2123— 005 observed with UVES-VLT clearly showing several 
H2 transitions, each comprising two distinct velocity features. Left: region 3097 - 3106 A with 
relatively low SNR. Right: region 3251 — 3263 A toward longer wavelength at higher SNR, also 
containing a line of HD. The upper trace shows residuals in the fit regions. The tick marks above 
the spectrum refer to the velocity components of the H2 and HD lines (green), modeled Lyman-a 
forest lines (blue), and Fe 11 lines (red). H2 spectral line identifications are shown at the bottom. 

llh observation in August and September 2008 at VLT delivered a spectrum reaching a 
signal-to- noise ratio of 65 per 1.5 km s~^ pixel at 3415 A, degrading to ~ 15 at the far 
UV end (3046 A). In this range many hydrogen lines are found, of which 90 H2 and 6 HD 
are useful to include in an analysis; two typical spectral regions are shown in Fig. [T| while 
the entire spectrum is reproduced in the additional material [H]. Individual spectra (14 
exposures of 2900s) were air-vacuum and heliocentrically corrected, rebinned and coadded. 

This high-quality VLT spectrum of J2123 was subjected to a comprehensive fitting pro- 
cedure with a direct minimization of A/i//i to the entire spectrum, in a similar manner as 
in [ini [12] . This procedure has many advantages over a line-by-line fitting procedure as used 
in other studies [71 [T31 [S]. The former allows more objective selection of spectral regions 
to be fitted and lines that are partly blended by Lyman-a forest lines can be modeled and 
explicitly included in the fitting routine. By this means the amount of information on H2 
absorption extracted from the spectrum is effectively increased; however, the regions where 
no H2 absorption occurs are cut out of the fitting routine to avoid contaminating the x^- In 
total, 68 regions cover the 96 H2/HD lines included in the fit. Additional regions containing 
some 39 H2/HD lines were excluded, because lines were too weak, or located within a satu- 
rated (Lyman-a) absorption region, or were contaminated with unidentified narrow spectral 
features, i.e metal lines. 

A generic model spectrum is produced from the accurately known transition wave- 
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lengths for H2 [8] and HD lines [9]. Intensities derive from multiplying calculated oscillator 
strengths [T3] with column densities N{J) for each rotational state as fitting parameters. 
This mask spectrum is then convolved threefold, with a Lorentzian for the damping coeffi- 
cient, Fj, for each line [16], a fixed width parameter for the instrumental profile of the UVES 
spectrograph (0.8" resulting in a resolving power of 53000), and a Gaussian fit parameter, 
b, for the Doppler width. In this way the H2 and HD lines in the spectrum are tied, leading 
to a reduction in the number of parameters to be adjusted in the fitting procedure. 

This spectrum is replicated to accommodate the velocity structure of the absorbing cloud, 
starting with the two features split by ~ 25 km/s visible in Fig. [!} Subsequently, additional 
velocity components (VCs) are added, thereby increasing the number of fit parameters by 6vc 
and Nvc{J) for all VCs. Only after a full II2/IID velocity structure is imposed, all Lyman- 
forest and metal lines are assigned and simulated properly, and the continuum and zero 
levels are adjusted for all local fitting regions, a possible variation of the proton-to-electron 
mass ratio A/i/yU is included as a final fitting parameter. The VPFIT algorithm [T7] uses 
Eq. ([2]) and the Ki coefficients |5l El dU] directly. In addition to Lyman-a forest lines, metal 
lines falling in the H2 regions are also co-fitted. As one example we mention the 18 velocity 
components of the Fe ll line, weakly overlapping the II2 L3R0 line at 3252 A [llj; its velocity 
structure is tied in the fit to the Fe ll line observed at 7952 A (Aiab = 2600 A). 
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FIG. 2: (Color online) Composite residual spectrum composed of the 25 cleanest H2 lines for fits to 
2, 3 and 5 velocity components. Upper panel: residuals. Lower panel: the spectrum as simulated 
from the final parameters resulting in the 2 and 5 VC models. The tick marks indicate the fitted 
positions in 2 and 5 VC models; in the 5 VC model two VCs lie at almost the same velocity and 
appear as a single tick mark on this plot. 
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FIG. 3: (Color online) Results from various fits, a) resulting value for Afi/fi for progressively 
increasing the number of velocity components; there are two results for 3VC fits, obtained by 
putting the 3'^'^ component in both distinct velocity features, which are horizontally offset for 
clarity; b) the Xu value; c) result on Aji/fi from fits to 7 subclasses of spectral lines. 

One of the advantages of fitting all lines simultaneously is that the velocity structure 
of the absorbing system can be robustly investigated. For example, a composite residual 
spectrum was constructed from the 25 cleanest absorption lines, all shifted to a common 
redshift {z = 2.059400), as shown in Fig. |2} It demonstrates that significant H2 absorption is 
left in a systematic way in the two-component fit; four and a five-component fits clearly yield 
better representations of the spectrum. Besides visual inspection of residuals, the procedure 
of adding velocity components can be quantified by calculating the value of per degree 
of freedom, x^, as is shown in Fig. [sj xl decreases until a possible 6**^ VC, which is then 
rejected on the basis of statistical grounds. Application of a modified version of the Akaike 
Information Criterion for finite sample sizes, yields a value of A^jc ~ 90 (for definition 
and evaluation see [IE]), thus establishing an objective ground for adding a 5*^^ VC, and for 
leaving out the 6*^ VC. The result of the 5-component fit is thus obtained as the fiducial 
statistical outcome with Afi/fi = (8.5 ± 3.6) x 10~^. 

The robustness of the statistical analysis was further investigated by verifying that the 
result was not produced by a limited subclass of lines. A fit was performed whereby a 
different value of Afi/fi was fitted simultaneously to each of the rotational states in hydrogen, 
Il2(J), and for HD. The results shown in Fig. ^ demonstrate the consistency for the value 
of Afi/fi in this procedure. In the model description the physical assumption was made 
that the Doppler widths, 6vc (for each velocity component z^c), were independent of the 
rotational state; if, however, a fit is performed with varying Doppler widths bv{J), a Ayu//i 
value is returned differing by only 5%. 
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The analysis strategy aimed to stick as close as possible to the known molecular physics 
of H2 and HD. The damping parameters, Fj, from (TB] describe the wings of the saturated 
lines (as in Fig. [T]) well. The oscillator strengths, /j, were kept fixed to the molecular physics 
values [15], thereby keeping the number of fit parameters as low as possible. However, in the 
region 3345 — 3415 A, covering the LO and LI Lyman lines, the calculated underestimate 
the line intensities by some 50%. Hence, in the fits the oscillator strengths in this range were 
adjusted, similarly as in [10], where this phenomenon was attributed to the O lv/Lyman-/3 
emission feature of the background quasar. 

As for systematic effects on the resulting value of Ayu/yU, adverse influences on the wave- 
length calibration are of particular importance. Wavelength calibration of the composite 
spectrum was derived from Th-Ar frames attached to each science exposure, where use was 
made of the improved Th-Ar line selection algorithm rejecting blended and weak lines [T9] . 
Wavelength calibration residuals of Th-Ar lines, used for calibrating the echelle grating are 
typically 70 m/s. However, since the wavelength scale of each echelle order is derived from 
typically more than 11 Th-Ar lines and many H2 transitions are fitted for each order, this 
effect largely averages out, and long-range variations of the wavelength scale are limited to 
30 m/s [19]. In view of a span in Ki of ~ 0.05, a shift of 30 m/s translates into a maximum 
variation in A/x/yU of 2.0 x 10"^; here velocity shifts are expressed as Av = cKiAfi/fi. 

In addition, intra-order distortions within individual echelle orders, first studied for the 
HIRES-Keck grating system [2S] and later for UVES-VLT [27], may result in peak-to-peak 
velocity shifts of up to Av = 200 m/s between transitions. This phenomenon was simulated 
by producing a counter-distorted spectrum by shifting the center of all echelle orders by 
— 100 m/s and imposing shifts up to -1-100 m/s at the order edges. From a fit to this 
spectrum it is estimated that this phenomenon produces a possible systematic effect of 0.7 x 
10^®. Systematic effects from spectral re-dispersion occurring in the re-binning process when 
averaging over the 14 science exposures are less than 0.2 x 10^^. Finally, drifts in temperature 
and atmospheric pressure between the quasar and Th-Ar calibration exposures are < 1 K and 
< 1 mbar, respectively, leading to a possible error of < 0.7 x 10^^. Taking the possible shifts 
in quadrature this results in an estimate of the systematic error of 5syst(A/i//i) = 2.2 x 10^^. 
Uncertainties in the Ki coefficients may lead to uncertainties in A^/fi at the ~ 0.01 x Afi/ n 
level [H], i.e. they are negligible. 

A A/i-effect can be mimicked by a longe-range distortion of the wavelength scale due 
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FIG. 4: (Color online) Current constraints on A/i//i on a cosmological time scale. Black points 
for z < 1 refer to NH3 values from objects B0218 [23121] and PKS1830 [25J; blue diamonds refer 
to Ref. [12]; green triangle to Ref. [10], i.e. the HIRES-Keck analysis of J2123. The grey circle 
refers to a recent re-analysis of Q0347 [H] . The red square is the present result of our UVES-VLT 
analysis of J2123. All uncertainties are at la. Independent constraints on the same object have 
been offset in redshift for clarity. 

to the correlation of i^-coefficients with wavelength [6J. Separate fitting of the range A < 
3230 A, where Werner and Lyman lines exhibiting strongly differing Ki values are found 
alongside each other, returns a value within 10% of the fiducial value, although with a 
larger uncertainty margin. This can be considered as an internal consistency check of the 
wavelength scale. 

The statistical and systematic analyses lead us to a final result for the UVES-VLT H2 
absorption spectrum toward J2123 (uncertainties representing la): 



Afi/fi = (8.5 ± 3.6stat ± 2.2,y,t) x 10" 



(3) 



This may be compared to the result from the analysis of the HIRES-Keck spectrum for the 
same object pUj : 

A/i//i = (5.6 ± 5.5stat ± 2.9syst) x 10"^ (4) 

These values are well within la of each other and may be averaged to yield Ayu/yU=(7.6± 
3.5tot) X 10~^. The agreement obtained between both investigations of the same high redshift 
H2 absorption spectrum from the two major large telescopes VLT and Keck enhances the 
confidence in the statistical and systematic uncertainty analyses. They are indicative of the 
error sources being suitably taken into account for both instruments. 
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Spectral observations of highly redshifted H2 absorbers are limited. Of the many thou- 
sands of known quasar systems, of which some 1000 have been classified as damped-Lyman-a 
(DLA) systems, in only some twenty is molecular hydrogen detected, but not under condi- 
tions allowing for a detailed analysis of /i-variation at a competitive level. So far only high 
quality spectra from H2 absorbing systems toward Q0347— 383 at Zabs = 3.02 and Q0405— 443 
at ;zabs = 2.59 [I2j and Q0528-250 at ;Zabs = 2.81 [I2j have been observed with the Very 
Large Telescope to yield a tight constraint on Afi/fi. Recently a VLT-reinvestigation of 
Q0347— 383 was also reported j5T]. In contrast, for detecting drifts in a, spectra of atomic 
species (Mg I, Mg ll, Fe ll, Zn ll, Cr ll. Si ll, etc.) in almost 300 absorption systems have 
been analyzed [20]. For H2, a spectrum of J1337+315 at z^hs = 3.17 was analyzed [22] . 
but due to the low column density (logA^(H2) = 14.1) a constraint on Afi/fi only at the 
non-competitive 10~^ level was found. 

In Fig. |4] the now existing information on high redshift absorbing H2 systems is collected; 
the tight constraints from the ammonia method [2SH2S] for z < 1 are included as well. The 
status result of a possible variation of the proton-to-electron mass ratio can be expressed 
as: there exists no firm evidence for a drifting /i on a cosmological time scale, which is 
constrained by Afi/fi < 1 x 10^^ at redshifts in the range z = 2 — 3. Quantitatively, 
averaging the six results for the four high-redshift H2 absorption systems yields A/i//i = 
(5.2 ± 2.2) X 10~^, which is a slight indication (at 2.3 a level) for a larger fi at high redshift, 
but further systems should be analyzed to warrant such a claim. As for the proportionality 
relation of Eq. ([T]), combined with the finding of Aa/a = (—5.4 ± 1.2) x 10^^ [20], some 
significant drift on fi exceeding the present constraint would be expected. The result may 
also be interpreted as the specific assumptions of Grand Unification underlying Eq. ([T]) to be 
unjustified. Indeed, while the majority of the GUT scenarios predict |-R^q| to be large, some 
special cases are reported yielding a smaller proportionality constant, i.e. \R^a\ < 5 [3]. 
The present results constraining Afi/fi may thus be used to discriminate between various 
possible GUT scenarios. 

This work is based on observations carried out at the European Southern Observatory 
(ESO) under program ID 81.A-0242 (PI Ubachs), with the UVES spectrograph installed 
at the Kueyen UT2 on Cerro Paranal, Chile. The authors acknowledge fruitful discussions 
with J. King and J.K. Webb (Sydney). MTM thanks the Australian Research Council, and 
WU thanks the Netherlands Foundation for Fundamental Research on Matter (FOM) for 
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